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ABSTRACT 

We present spin-resolved spectroscopy of the accreting white dwarf binary V455 And. With 
a suggested spin period of only 67s, it has one of the fastest spinning white dwarfs known. To 
study the spectral variability on the spin period of the white dwarf, we observed V455 And 
with 2 s integration times, which is significantly shorter than the spin rate of the white dwarf. 
To achieve this cadence, we used the blue arm of the ISIS spectrograph at the 4.2-m William 
Herschel Telescope, equipped with an electron multiplying CCD (EMCCD). Strong coherent 
signals were detected in our time series, which lead to a robust determination of the spin period 
of the white dwarf (Pspin = 67.619 ± 0.002s). Folding the spectra on the white dwarf spin 
period uncovered very complex emission line variations in Hy, He I X4412 and He II X46S6. 
We attribute the observed spin phase dependence of the emission line shape to the presence 
of magnetically controlled accretion onto the white dwarf via accretion curtains, consistent 
with an intermediate polar type system. We are, however, not aware of any specific model that 
can quantitatively explain the complex velocity variations we detect in our observations. The 
orbital variations in the spectral lines indicate that the accretion disc of V455 And is rather 
structureless, contrary to the disc of the prototype of the intermediate polars, DQ Her. This 
work demonstrates the potential of electron multiplying CCDs to observe faint targets at high 
cadence, as readout noise would make such a study impossible with conventional CCDs. 

Key words: binaries: close - stars: individual (V455 And) - novae, cataclysmic variables - 
white dwarfs - stars: magnetic fields. 



1 INTRODUCTION 

V455 And, also known as HS 2331+3905, is a grazingly eclips- 
ing short period cataclysmic variable (Porb = 81.08mins) identi- 
fied from the Hamburg Quasar Survey whose vari ability is very 
complex (Araujo-Betancor et al. '2005*: 'Pv rzaslEoill) . It contains a 
rapidly s pinning white dwarf accreting from a low mass compan- 
ion star. lArauio-Betancor et al ] (120051) found two closely spaced 
signals at 67.2 s and 67.6 s, and labelled the strongest of the two, 
at 67.2 s, as the white dwarf spin period. By studying the stabil- 
ity of bot h signals in a longer photometric dataset spanning 25 d, 
iGansickd (12007) concluded that the signal at 67.6 s must be related 
to the true spin period. The shorter period signal is more likely 
to be a beat between the spin period and the ^ 3.5 h periodicity 
detected in spectroscopy at some epochs jArauio-Betancor et al.i 
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I2OO5I : iTovmassian et alll2007l) . V455 And also shows variability 
which could be interpreted as white dwarf pulsatio ns. The sys- 
tem e xhibited its first known super-outburst in 2007 (iMatsui et al.l 
2009). The temperature of the wh i te dw arf is 10500 ± 500K in 
quiescence dArauio-Betancor et al ] E005D . Dwarf nova outbursts 
heat up the white dwarf, which make V455 And an important 
test bed for our understanding of the white dwarf instability strip 
(Szkody et al. 2010; Silvestri et al. 2012). 

In this study, we focus on fast spectroscopy resolving the spin 
cycle of the accreting white dwarf. A strong modulation on the 
white dwarf spin is expected when the white dwarf is sufficiently 
magnetic to control the accretion flow onto its surface, funnelling 
the material onto its poles. Cataclysmic variables with a magnetic 
white dwarf whose magnetic field is not strong enough to synchro- 
nise the white dwarf spin with the orbit of the system are known 
as intermediate polars (IPs). In most of the known IPs, including 
V455 And, the region where the accretion flow is controlled by 
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the magnetic field is limited in radial extent from the white dwarf 
and the bulk of the accretion flow still proceeds via a traditional 
accretion disc that is steadily fed by the mass donor star. Studies 
of the spectral line variations of IPs on the spin period have been 
performed for a se ries of IPs with slowly spinning WDs (see e .g., 
iHelHer et al.lll987l for a spin trail of Ex Hya, Hellier e taDll99Ql for 
FO Aqr, Helher et al. 1991 for AO Psc, Helher 1997 for BG CMi 
and PQ Gem, and Helher 1999 for V2400 Oph and V1025 Cen). Of 
the IPs with faster spinning WDs, however, only spin var iations of 
AE Aqr dWelsh et al."l998': 'Reinsch & Beuermannl ll994l) and DQ 
Her teloemen et al. 2010; Martell et al. 1995) have been studied. 
The observed changes in the DQ Her spectra remain largely unex- 
plained. No significant spectral line variations have been detected 
on the spin period of AE Aqr. This could be e xplained by the fact 
that the system is probably a propeller system teracleous & Horna 
[l996; Wynn et al. 1997) in which most of the material lost from 
the companion gets expelled from the system rather than being ac- 
creted on the white dwarf. 

The earlier papers on IPs f see |Pattersonlll994l for a review), 
suggested variations in the light curve on the spin period to arise 
from X-ray reprocessing in the disc or bright spot. Reprocessing 
of X-rays from the magnetic poles of the WD would lead to sinu- 
soidally shaped brightness enhancements in the spin trails, as the 
beam sweaps around and illuminates different parts of the disc at 
different spin phases. More recently, it has become clear that mag- 
netically controlled accretion via accretion curtains is likely to have 
an important influence on the observed light, by its changing ori- 
entation (and hence visibility) during a spin cycle, and possibly by 
acting as an extra reprocessing region. The idea of curtain- like ac- 
cretion flows from the disc to the magne t ic pole s of the white dwarf 
was originally proposed by Rosen et alj (Il988h . 

To better understand the accretion geometry near rapidly spin- 
ning white dwarfs, we targeted V455 And which, given its ex- 
pected spin period of only ~ 67s, is one of the most rapidly 
spinning white dwarfs know n after AE Aqr (Pspin = 33.1s, 
Patterson|[T979l:lde Jager eTaD[T994 ) and V842 Cen (P^pin = 56.8 s. 



at high cadence. They are therefore also known as Low Light Level 
CCDs ( L3CCDs). QUCAM 2 was used earlier on the red arm of 
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Woudt et al.ll2009l) . The required observations are technically chal- 



lenging given the need to obtain relatively high spectral resolution 
spectroscopy with very short exposure times. We describe the in- 
strumental setup involving an electron multiplying CCD in Section 
[21 detail our custom spectral extraction procedure in Section (3] fol- 
lowed by our data analysis and interpretation. 



ISIS bv lTulloch et al] (l2009l). M ore information on QUCAM2 can 



be found in iTuUoch & DhillonI (|201Q|) . We paired QUCAM2 with 
the R1200B grating and a 1" wide longslit to cover 4300 to 4700 A 
at a spectral resolution of 1.5 A. In total 15 796 blue arm spectra 
were taken with an integration time of 2.0 s. The observations are 
summarized in Table [T] The slit was held at a fixed sky- angle to 
permit the simultaneous observation of a brighter comparison star, 
allowing us to monitor slit-losses and telescope tracking. 



3 DATA REDUCTION 

STARLINj|3 routines and the MOLLY software packages written by 
TRM were used for the data reduction. 



3.1 Gain determination 

The gain of an EMCCD is highly dependent upon the voltage that is 
applied to the serial register that amplifies the signal before readout. 
We determined the gain of the EMCCD from the overscan regions 
of the science spectra. Since these have no illumination, the number 
of electrons is zero in most pixels except in a few pixels because of 
clock induced charges. The probability distribution for the output 
of 1 input electron is approximately 



P{x)n= 



exp(-x/g) 



IBasden et al.ll2003h . The histogram in Fig. [T] shows the logarithm 
of the distribution of the pixel values of the overscan regions of all 
science frames. The Gaussian distribution around results from the 
pixels with zero charge. The linear part results from the pixels with 
a single electron charge. The number of pixels with a higher charge 
is negligible. The slope of the linear part is — 1/g, from which we 
derive g = 170 ± 8 for the combined data of the three nights. There 
was no evidence for significant gain variations between different 
nights, thus we use the average gain for all extractions. This value 
is slightly higher than the nominal gain of 160 mentioned in the 
instrument manual, but longer-term variations at this level are to be 
expected. 



2 OBSERVATIONS 

V455 And was observed with the double armed ISIS spectrograph 
mounted on the 4.2-m William Herschel Telescope on La Palma 
(Spain), on the nights of 11, 12 and 13 July 2008. Since the inte- 
gration time can only be of the order of a few seconds, the signal 
of V455 And would be buried under the readout noise of a conven- 
tional CCD. The camera we used, QUCAM2, is equipped with a 
Ik X Ik Electron-Multiplying CCD (EMCCD) in which the signal 
is amplified before readout, such that even the signal of a single 
photon dwarfs the readout noise. In addition, the CCD has a frame 
transfer buffer which allows on chip storage of an image such that 
an exposure can be started while the previous is still being read 
out. This makes it possible to observe time series of spectra with 
a negligible dead time of about 12ms between two frames. EMC- 
CD s are more efficient than conventional CCDs when the number 
of cou nts per pixel is smaller than the readout noise squared (iMarshI 
l2008h . which makes them particulary useful to observe faint targets 



3.2 Extraction and calibration of curved spectra 

The spectra of V455 And were slightly curved and tilted. There- 
fore, the spectra were extracted along a 3rd order polynomial fit 
to the spectrum (trace) following .Marsh (.1989.) . Because the ad- 
vantage of EMCCDs is largest when working at low count levels 
(Marsh 2008), the science spectra typically have low count levels. 
This prevents one from accurately fitting the trace of the curved 
spectra. If the position of the spectra on the CCD is constant in 
time, this problem can be overcome by stacking lots of frames. 

Unfortunately, the positions of the object spectra were found 
to be wobbling along the spatial axis on timescales of roughly a 
minute with a semi- amplitude of 2 to 3 pixels, equivalent to about 
0.5 arcsec. To deal with this extra complication, we implemented a 

^ STARLINK is Open source software and can be obtained from 
http :// starlink.j ach . hawaii . edu/ starlink . 

^ MOLLY and DOPPLER are available for download at 
http ://deneb. astro .Warwick, ac.uk/phsaap/software/ . 
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Table 1. Overview of our spin-resolved spectroscopic observations of V455 And, performed with an EMCCD at the WiUiam Herschel Telescope. 



Date 


UT 




Conditions 


Seeing 


Instrument 


Grating 


Exp. time 


Number 


11-12/07/2008 


02:02 


- 03:59 


clear 


0.8 to 1.8 arcsec 


QUCAM2 at ISIS blue arm 


R1200B 


2.0s 


3102 


12-13/07/2008 


01:13 


- 05:39 


clear 


^ 0.7 arcsec 






2.0s 


7020 


13-14/07/2008 


01:46 


- 05:25 


clear 


0.5 to 0.8 arcsec 






2.0s 


5674 




200 300 400 

ADU - bias 

Figure 1. Histogram with 400 bins of the pixels values around the bias level 
of the pixels in the overscan region of the science spectra. All pixels are 
expected to have had a charge of or 1 electron. The Gaussian distribution 
centred at represents the pixels with a zero charge. The pixels with a 
charge of one electron give rise to the tail, that extends to very hight count 
values. The slope of the linear fit to the tail is — 1/g. 



multistep re duction strategy whi ch replaces the standard reduction 
techniques aTulloch et alJ (l2009h used to reduce QUCAM2 data. 
First, we determined a 3rd order polynomial trace and a profile for 
optimal extraction for both stars from a stacked image. Then, we 
tweaked the first order parameter of the trace of the (brighter) com- 
parison star on each frame. Subsequently, we shifted the frames to 
align the spectra of the comparison star, and thus also the spectra 
of the target assuming a constant offset between the two on short 
timescales. From the shifted frames, a new stack was made which 
was used to recompute the traces and profiles. These traces and 
profiles were used to extract the spectra of the two stars from the 
original frames, tweaking the first order parameter of the trace of 
the comparison star and applying the same shift to the trace of the 
target star. The sky regions that were used for the sky subtraction 
were shifted by the same number of pixels. 

Arc spectra were taken every ~ 2 hours and wavelength cali- 
brated by a third order polynomial fit to the lines achieving an RMS 
of ~ 0.05 A. The wavelength scale of the V455 And spectra was 
determined by linear interpolation between bracketing arc spectra. 
Spectra of HZ44 were taken every night and used to flux-calibrate 
the V455 And spectra. Slit losses of individual exposures were ac- 
counted for by comparison of the flux level of the spectrum of the 
in- slit comparison star with a wide slit spectrum. 



4 ANALYSIS 

Fig. [3 shows an example single exposure (grey) and the average 
(black) spectrum of V455 And. Double peaked line profiles of Hy 
M340 and He I X447\, characteristic for accretion disc emission, 
can clearly be seen. He II X46S6 is visible as well, but is much 
weaker. 




Wavelength (A) 

Figure 2. The average (black) blue arm spectrum of V455 And overplotted 
on a single spectrum (grey). 



4.1 Frequency analysis 

To search for oscillations in the emission lines, we present a Lomb- 
Scargle periodogram (Scargle 1982; Lomb 1976) of the net line flux 
in Hy in Fig. [3] The net line flux was obtained by integrating the flux 
in the line (above the continuum level) over a region of 60 A cen- 
tered at Hy. Our short exposure times allow us to resolve line flux 
modulations out to several thousand cycles/day. The highest peaks 
are found at the orbital frequency of 17.5 cd~^, and its one-day 
aliases (the window function is shown in grey in the insets, cen- 
tered at the highest peak). In addition, the highest amplitude peak 
in the second inset is at 1277.8 cd~^. This coherent signal would 
be most naturally identified as the spin period of the white dwarf, 
as IPs often show a spin-locked signal in their emission lines. Fur- 
thermore, no peak is found ab ove noise level at 1284.7 c d~^ , which 
lArauio-Betancor et al ](|2005|) suggested to be WD spin frequency 
but is more likely a beat frequency, as explained in Section [T| The 
third inset zooms in on the region around 2555.6 cd~^ the first 
harmonic of the spin period whose peak power is higher than the 
fundamental. The variation of the spectra on the orbital period is 
further discussed in Section |4l2l and the variability on the spin pe- 
riod in Section |431 

We do see excess power in the emission line power spectrum 
near the nominal pulsation period range (~ 300c d~^). This implies 
that the pulsation source is not a pure continuum source, but must 
have a variability component in the Hy line. This could be caused 
by variability in the underlying B aimer absorption profile from a 
pulsating white dwarf, but is also consistent with a variable contri- 
bution from the line emission sources. Our data cannot distinguish 
between these two given that both the white dwarf absorption line 
and the accretion powered emission line have comparable widths 
and thus, at all velocities where w e detect variability, overlap (see 
Fig. 12 in lArauio-Betancor et al.ll2005h . 

We have also checked the variability in the continuum part 
of the spectra (between 4520 and 4670 A). We found variability at 
the same frequencies as in the Hy line, but at lower amplitude, ex- 
cept for the variability at 1277.8 c d~^ which is not detectable above 
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Frequency (c/d) 



Figure 3. Scargle periodogram of the flux in Hy above the continuum. In 
the insets, the window function is plotted in grey, shifted to the frequency 
of the strongest peak. The highest peak in the left inset is 18.5 c d~^ which 
is the one day alias of the orbital frequency peak at 17.5 cd~^ which is 
almost equally strong. In the second inset, the highest amplitude is found 
at 1277.8 cd~^ which is identified as the spin period of the white dwarf. 
The peak in the third inset, at 2555.6 c d~\ is the first harmonic of this spin 
period. 

the noise level. This is in agreement with the photometric studies 
mentioned earlier, which find a much lower variability amplitude at 
1277.8 c d~hhan at the first harmonic of this frequency. 

4.2 Variations on the orbital period 

We rebinned the spectra on a constant velocity scale of 28kms~^ 
per pixel and folded the spectra on the orbital ephemeris of 
lArauio-Betancor e t~ar (2005). The resulting orbital trails for Hy, 
He I XAA12 and He II M686 are shown in the upper panels of Fig.|4] 
Apart from the double peaked profile emitted by the accretion disc, 
a strong bright spot contribution is seen as an S-wave in the trails . 
A weak WZ Sge-like 'bright spot shadow' dSpruit & Ruttenlll998h 
can be seen in the red-shifted part of the lines at orbital phase 
^ 0.25. Around orbital phase 1, a shallow eclipse is observed. In all 
three lines, the red-shifted peak appears to be brighter than the blue- 
shifted peak at all orbital phases. This difference is probably due to 
a slightly warped or slowly precessing accretion disc of which the 
red-shifted surface was better visible and/or brighter at the time 
of our observations compared to the blue-shifted side. Comparing 
night-by-night averages reveals no significant evolution in this line 
asymmetry across our three nights. 

The bottom panels of Fig. |4l show the Doppler maps of the 
tra ils, produced using the maximum entropy method as presented 
bv lMarsh&Homel(ll988h and implemented in the DOPPLER^ pack- 
age. The three maps are not significantly different. The accre- 
tion disc emission maps onto a diffuse ring. The disc appears 
to be structureless, as is the case for most IPs but contrary to 
the prototype DQ Her, which wa s found to have spiral den- 
sity structures (iBloemen et alJEoiol) . The dot around (Vx,V7) ^ 
(—900, 100)kms~^ represents the strong brightspot emission com- 
ponent which is hot enough to contribute significantly to Hell. No 
line emission from either the white dwarf or the secondary star is 
detected. These maps show that the bulk of the line emission orig- 
inates from an extended accretion disc, implying that any magneti- 
cally controlled flow is confined to the zone near the white dwarf. 

4.3 Variations on the WD spin period 

To exploit our EMCCD setup and study the changes in the line 
profiles at the white dwarf's spin period, the spectra were phase 



folded into 50 bins at the spin period of 67.619 s (1277.75 cd~^). 
An average spectrum was subtracted to highlight the spin modu- 
lation, similar to the procedure outlined in Bloemei Tet al.l feOlOh . 
The resulting spin trails are shown in Fig. |5]for the three spectral 
lines. The detection of a clear pattern i s strong evidence for the fact 
that, as concluded bv lGansickg (l2007h . 67.619 s is the true spin pe- 
riod because folding the spectra on slightly different per iods, such 
a s the 67.24 s period that was suggested by Arauio-Be tancor et alj 
i2005h . washes out the structure. In DQ Her there is still a debate 
whether the spin signal is at the spin period or its harmonic. Here 
we can confidently conclude that the folded trails robustly desig- 
nate 67.619 ± 0.002 s to be the rotation period of the white dwarf, 
with the power detected at the 1st harmonic caused by the two-sided 
pattern in Fig.|5] The uncertainty on the spin period was determined 
by comparing how well spin trails of the different nights phase up, 
after folding the data on various different periods that are slightly 
different from the optimal spin period value. 

The observed pattern is puzzling. Previously published spin 
trails of IPs with relatively slowly spinning IPs (see references 
in Section [T]) appear to be totally different to the trails we find 
for V455 And. Also the recently published spin trails of DQ Her 
(see Bloe men et al.ll20ig) . although also being an IP with a fast- 
spinning WD, are fundamentally different. Firstly, DQ Her's vari- 
ations in He II M686 are found to be much stronger in the red- 
shifted parts of the line than in the blue- shifted, which is clearly 
not the case for V455 And. The spin modulations can be detected 
out to large velocities and their amplitudes peak at velocities sig- 
nificantly beyond the disc peak. To illustrate this, we plot several 
representative normalised spin profiles in Fig. |6] These illustrate 
the amplitude of the observed spin signal as a function of velocity, 
with an average emission line profile provided as a reference. The 
structure in the spin trials and profiles cannot be explained by the 
'old school' hypothesis of X-ray reprocessing in the disc and the 
brightspot. If reprocessing of X-rays by the bright spot contributed 
significantly, we would expect a higher spin modulation amplitude 
in the spectra lines at the radial velocity of the bright spot, which 
would cause the spin pattern to be orbital phase dependent. Spin 
trails were therefore also produced using subsets of the spectra, 
taken at selected orbital phase intervals, but no variability of the 
spin trails over the orbital period could be detected. The complex 
pattern more likely results from magnetically controlled accretion 
near the white dwarf. Accretion curtains flowing towards the white 
dwarf would involve large velocities, switching from maximum to 
minimal radial velocity as they rotate along the line of sight across 
1/4 spin cycle. However, we do not know which exact geometry 
can explain the observed modulations. A thorough modelling ef- 
fort of the accretion curtain geometry and its optical properties will 
be required to quantitatively test this idea and see whether such a 
geometry could reproduce the observed spin modulations as a func- 
tion of radial velocity and spin phase. 



5 SUMMARY 

We observed the intermediate polar V455 And with the electron 
multiplying CCD 'QUCAM2' installed at the ISIS spectrograph at 
the William Herschel Telescope. We demonstrate the potential of 
this detector to observe faint targets at high cadence. With a con- 
ventional detector, spin-resolved spectroscopy of the IPs with the 
fastest spinning WDs would be hard if not impossible to achieve as 
readout noise would swamp any signal. We developed a strategy to 
reduce the faint spectra, making full use of a brighter in- slit com- 
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Figure 4. Orbital trails and Doppler maps for Hy, He I A,4472, He II A,4686. White is minimum flux, black maximum flux. A grazing eclipse can be seen at 
orbital phase 1. The three lines result in comparable Doppler maps that show emission of the accretion disc (ring) and the bright spot where the mass stream 
from the companion hits the disc (brighter dot at the ring). 
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Figure 5. Trails of normalised and mean subtracted line profiles, folded on V455 And's spin period of 67.619 s. The color scales are ± 0.5, 0.2 and 0.1 mJy 
compared to average. White indicates lower than average fluxes, black higher than average. 



parison star and correcting for instrument and telescope flexure. 
We were able to detect strong coherent signals in our time- series 
allowing us to robustly identify the spin period of the white dwarf 
to be 67.619 lb 0.0 02 s, which confirms the spin period reported 
bv lGansick3 (l2007 l). Furthermore, by folding our 15 796 spectra on 
this spin period, a complex emission line variation can be recovered 



resolving the spin modulation as a function of both radial velocity 
as well as spin phase. The observed variations are totally differ- 
ent from the results of previous observations of other intermediate 
polar s, including the canonical IP DQ Her. We believe that the ob- 
served patterns are evidence of magnetically controlled accretion 
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Figure 6. Five spin profiles for the Hy line normalised to show the spin 
modulation in percent compared to the average line profile at the top. The 
spin modulation peaks at velocities well beyond the disc peaks and can be 
detected out to large radial velocities. 



curtains near the white dwarf, but are not aware of any specific 
model that can reproduce our observations in detail. 
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